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ABSTRACT 

We present a model for the dispersal of protoplanetary disks by winds from either the central star 
or the inner disk. These winds obliquely strike the flaring disk surface and strip away disk material by 
entraining it in an outward radial-moving flow at the wind-disk interface which lies several disk scale 
heights above the mid-plane. The disk dispersal time depends on the entrainment velocity, = ec^, 
at which disk material flows into this turbulent shear layer interface, where e is a scale factor and 
Cs is the local sound speed in the disk surface just below the entrainment layer. If e 0.1, a likely 
upper limit, the dispersal time at 1 AU is ~ 6 Myr for a disk with a surface density of 10^ g cm~^, a 
solar mass central star, and a wind with an outflow rate = 10~* Mq yr~^ and terminal velocity 
Vyj — 200 km s^ . When compared to photoevaporation and viscous evolution, wind stripping can be 
a dominant mechanism only for the combination of low accretion rates (< 10~* Mq yr"-'^) and wind 
outflow rates approaching these accretion rates. This case is unusual since generally outflow rates are 
< 0.1 of accretion rates. 

Subject headings: planetary systems: protoplanetary disks — stars: winds, outflows — hydrodynamics 
— accretion, accretion disks 



1. INTRODUCTION 

A crucial timescale in planet formation is the timescale required for the loss of protoplanetary gas material which 
initially dominates the total disk mass. The outcome for a particular planetary system might be very different if the 
parent disk is dispersed faster or slower than in our solar system. However, relatively little attention has been paid to 
the processes responsible for dispersing protoplanetary disks. HoUenbach et al. (1994) set a theoretical framework and 
Shu et al. (1993) applied it to propose photoevaporation of the solar nebula as the disk dispersal mechanism capable of 
explaining the differences in envelope masses between the gas-rich giants, Jupiter and Saturn, and the gas-poor giants, 
Uranus and Neptune. HoUenbach et al. (2000) generalized the discussion, describing the variety of possible disk 
dispersal mechanisms. The dominant disk dispersal mechanism at the inner parts of the disk is viscous accretion onto 
the central star. However, this process becomes inefficient with time as the outer disk continuously expands to conserve 
angular momentum and the accretion rate decreases. Other possible disk dispersal mechanisms are planet formation, 
stellar encounters, stellar winds or disk winds, and photoevaporation by energetic photons. Photoevaporation is a 
process in which the surface of the disk is heated by stellar photons, resulting in a hydrodynamical flow, a slow wind 
(< 10 km s"^), back to the interstellar medium. HoUenbach et al. (2000) concluded that planet formation is a minor 
disk dispersal mechanism, and that the dominant mechanisms for a wide range of disk sizes are viscous accretion in 
the inner disk and photoevaporation in the outer disk. Therefore, planet formation must compete with these more 
efficient dispersal mechanisms. 

Handbury & Williams (1976) found that a stellar wind could not have removed the solar nebula. They argue that 
the stellar wind could only push the nebula to a finite distance determined by force balance and angular momentum 
conservation. However, they assumed that the stellar wind pushes the nebula as a whole and that the nebula maintains 
Keplerian rotation. Yun et al. (2007) studied the geometrical and thermal structure of the wind-disk interface for the 
specific case of a passive disk. However, they did not consider disk dispersal by wind stripping. Cameron (1973) 
suggested that the interaction of an outflowing stellar wind with the solar nebula would lead to significant mass loss, 
and Horedt (1978) and Elmegreen (1978) proposed models which did consider disk dispersal. Although both of these 
models predict significant mass loss, the manner in which the disk is dispersed is remarkably different. Horedt (1978) 
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Fig. 1. — (a) Mass flow from the wind (with velocity and incidence angle 7) and the underlying disk (with velocity ecs perpendicular 
to the mixing layer surface) to the wind-disk mixing layer. There is mass transport through the mixing layer in the direction of the local 
tangent angle, /3, at velocity Vm- (b) A mixing layer annulus of length 5s between R and R + 5R receives mass, momentum, and angular 
momentum from the wind (with velocity and incidence angle 7), the disk (with velocity eCs perpendicular to the mixing layer surface), 
and the inner annulus (with velocity v,n and incidence angle S/3), while it delivers mass and momentum to the outer annulus (with velocity 
). The local tangent of the mixing layer surface is given by dz/dR = tan/3. 

predicts that the wind drives an outward flow, while Elmegreen (1978) predicts, somewhat paradoxically, that the 
wind ultimately drives an inward flow. The difference between these two models is discussed in Elmegreen (1979). 
In particular, the wind-disk interface is given a priori in Horedt (1978), while it is calculated in Elmegreen (1978) 
by considering normal pressure balance. We follow the latter course here. Elmegreen (1978) showed that if the 
momentum is deposited where the wind strikes the disk, the addition of low angular momentum wind material causes 
the underlying layers to spiral inwards. As we will show below, if Keplerian rotation is assumed in the mixing layer, 
the net radial flow is inward. In this case, the wind causes disk dispersal by accelerating accretion onto the central 
star, as shown by Elmegreen (1978). 

We consider the alternative case of a shear mixing layer whose rotation is non- Keplerian but conserves momentum of 
wind and entrained disk surface material. The velocity shear is large and the tenuous surface material moves outward 
at speeds greater than the gravitational escape speeds. The disk material is expected to be entrained into the wind-disk 
interface in this case, and is carried outwards to the interstellar medium. The entrainment layer, or shear layer, carries 
with it a mixture of shocked wind material as well as entrained disk material. 

The observationally inferred wind mass loss rate in T Tauri stars can be as high as ^ 10~^ Mq yr~'^ (e.g. Cabrit et al. 
1990; Hartigan et al. 1995; White & Hillenbrand 2004). For comparison, the present mass loss rate of the Sun's 
stellar wind is 10~^^ Mq yr~^. The strong magnetic activity in young stars is capable of driving stellar winds 
with outflow rates < 10~^ M0 yr"-'^ (Decampli 1981). In young stars with accretion disks, the interaction of the 
rotating magnetic fleld with the accreting disk can generate outflows with mass loss rates < 10~^ Mq yr^^ for 
correspondingly high accretion rates < 10~® yr~^ (see review by Koenigl & Ruden 1993). It remains unclear 
whether the outflows are launched near the magnetospheric truncation radius (Shu et al. 2000) or over a wider range 
in disk radii (Konigl & Pudritz 2000). Our emphasis here is on the subsequent interaction of the wind, however 
launched, with the protoplanetary disk. Thus, hereafter we will simply refer to the stellar wind or the disk outflow as 
the central wind. 

The rest of the paper is organized as follows. §2 describes the theory, which incorporates pressure balance; and mass, 
momentum, and angular momentum conservation. §3 presents the results of applying the theory to protoplanetary 
disks. Finally, §4 summarizes the main results and discusses some of their consequences. 

2. MODEL 
2.1. Wind-disk mixing layer 

Following the analysis of HoUenbach et al. (2000), we assume that disk dispersal by the stellar wind occurs in a thin 
mixing layer at height z{R), where R is the disk radius. Disk material flows into this mixing layer with speed ec^, where 
Cg is the sound speed and e is the entrainment efficiency, as shown in Figure 1. This type of model has been applied to 
protostellar jets entraining circumstellar material by Canto & Raga (1991). Since the disk cannot react to changes at 
velocities faster than the sound speed, e = 1 is a natural upper limit to the possible value of the entrainment efficiency. 
Experimental results at high Mach numbers in the regime that is relevant for central winds suggest that e ~ 0.1 is a 
more likely upper limit (Canto & Raga 1991). We will consider entrainment efficiencies in the range 0.01 — 1. 

We denote the wind density, velocity, and incidence angle as pw, Vw, and 7 respectively; the disk density and 
Keplerian angular velocity as pd and f2fc; and the mixing layer surface density, velocity, inclination, and angular 
velocity as Em, Vm, P, and fim respectively. The wind incidence angle, 7, and the mixing layer inclination, /?, are 
related by 

a ^1 + 0, (1) 

where tanO = z/R (Figure 1). 

In Appendix A, we derive equations for pressure balance, and mass, momentum, and angular momentum conservation 
in the mixing layer. For the benefit of the reader, we repeat these equations here: 
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Fig. 2. — Relationship between the curvature of the mixing layer and normal pressure balance given by equation (2). The wind pressure 
is given by = Pwf^ sin^ 7 and the total disk pressure is given by Pa = pdCs(l + t'^). 
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d cos 13 d 

— (Emf^m) = Pd<^Cs^k {R^^mVm^m) ■ (5) 

Equation (2) describes normal pressure balance in the mixing layer. On the right-hand-side (RHS) of this equation, 
the first and second terms correspond to the normal pressures of the wind and the disk respectively, while the third term 
describes the centrifugal force associated with the mixing layer curvature. A similar, so-called centrifugal correction, 
term has been considered in previous studies (Hayes & Probstein 1966, p. 137, eq. 3.2.7; Canto 1980, eq. 10; 
Hartmann & Raymond 1989, eq. 7b; Wilkin & Stahler 1998, eq. 34). While these studies assume that the tangential 
momentum of the mixing layer is given by the total accumulated tangential momentum of the wind, we calculate it 
self-consistently by solving for E^, Vm, and /3 from equations (2)-(5). 

Figure 2 illustrates the relationship between the mixing layer curvature and normal pressure balance given by equation 
(2). If the wind pressure, p^v'^sm'^^, is equal to the total disk pressure, Pdcl{l + e^), the curvature, cos (3d[3 / dR, is 
zero, as shown in Figure 2a. If the total disk pressure is larger than the wind pressure, the curvature is positive and 
the mixing layer becomes concave on the side of the wind (i.e. the mixing layer surface flares outward with increasing 
radius), as shown in Figure 2b. Conversely, if the wind pressure is larger than the total disk pressure, the curvature 
is negative and the mixing layer becomes concave on the side of the disk, as shown in Figure 2c. 

Equation (3) describes mass conservation. At a given radius, disk material mixes not only with the wind but also 
with the wind material already mixed with the disk gas at smaller radii. On the RHS of this equation, the first and 
second terms describe mass input from the wind and the disk respectively, and the third term describes mass transport 
through the mixing layer. If we ignore the wind and disk mass input, equation (3) is similar to the standard continuity 
equation in cylindrical coordinates, 

(!:„.«„..) = 0. (6) 

The factor of cos [3 in our continuity equation (3) arises because of the non-zero, variable, inclination of the mixing 
layer. Note that if we ignore the wind and disk contributions, equations (3) and (6) are equivalent for /? = 0, as 
expected. 

Equation (4) describes tangential momentum conservation. The first and second terms on the RHS of this equation 
correspond to the wind contribution and the tangential momentum transport through the mixing layer respectively. 
Once again, the factor of cos/3 arises because of the non-zero, variable, inclination of the mixing layer. 

Equation (5) describes angular momentum conservation. The first term on the RHS of this equation corresponds to 
the angular momentum contribution from the underlying disk, while the second term describes the angular momentum 
flux through the mixing layer. We assume that the central wind does not carry appreciable angular momentum. 

Combining equations (3) and (5) yields 

_ [Pd<^Cs (r^fc - - T.mdt^m " ^mPwVw sin 7] 

"""" i]™cos/3afl(a„i?2) • ^ > 

If we assume Keplerian rotation for the mixing layer (ri,„ = fife and dt^m — 0); equation (7) shows that the mass 
flow must be toward the central star {vm < 0), which is essentially in agreement with the result of Elmegreen (1978). 
As we will show in the next section, outward flow is possible if we relax the assumption of Keplerian rotation in the 
shearing mixing layer at the interface of wind and disk. 
We can write the disk mass loss rate per unit time and disk gas surface area as 

^.s^'-^, (8) 
cos/3 



where the factor of two accounts for mass loss from both sides of the disk and the factor of cos /3 accounts for the 
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inclination of the mixing layer. The characteristic dispersal time at a given radius due to mass loss to the mixing layer 
can be written as 

2.2. Disk 

We assume hydrostatic equilibrium in the vertical direction to write the disk density as 

^ S],_exp [-zV(2g^)] 
^'^^ ' '''^ y27riJerf[z„/V2i/] 

for z < Zm, where is the gas surface density of the disk, z,„ is the mixing layer height, erf(a;) is the error function, 
and H is the disk scale height. The disk scale height is given by 

where kg is Boltzmann's constant, is the disk temperature at radius R, fi — 2.34 to/^ is the disk mean particle 
mass, and rriH is the mass of a hydrogen atom. For a thin flat disk, cx i?^'^/'*, while for a disk that flares outward 
(Eq. 11), a larger fraction of the stellar flux is intercepted and Td oc R^^^^ (Kenyon & Hartmann 1987). We assume 
a midplane disk temperature power law distribution, — 100 K{R/AU)^^^^ (D'Alessio et al. 1998). 

In the so-called "a-disk" theory (Shakura & Sunyaev 1973), the viscosity, v, is parameterized as v = aCgH, where 
a is a dimensionless constant parameter. For our adopted disk temperature profile, viscous diffusion causes the 
surface density profile to approach S^j oc R~^, regardless of specific initial conditions (Lynden-Bell & Pringle 1974; 
Hartmann et al. 1998). We adopt a power-law for the radial dependence of the disk gas surface density, 

E,(i?) = 10^ g cm-2 \ (12) 



The characteristic viscous evolution time scale is given by 

R^ ( R 

^ — ^0.4 Myrf — 



Comparison of predictions from similarity solutions with observed accretion rates and disk sizes suggests a ~ 10~^ 
(Hartmann et al. 1998), while the observed semimajor axis distribution of extrasolar planets suggests a ~ lO""' 
(Ida & Lin 2005). We will consider a values in this range. 
The accretion rate can be written as 

A/_ = 37r^.S, ^ 3 X 10-« Mq yr-^ f , ^° \ ( ( , (14) 

where Sq and Tq are the surface density and disk temperature at 1 AU. The radial velocity of the accretion fiow is 2 
cm s" -^a/lO"^)- Hence our assumption that only material and not radial momentum crosses into the mixing layer. 
The inner boundary is chosen such that the wind velocity is equal to the escape velocity from the star at this radius: 

= ^ ^ 0.2 AU (^) ( ,) " . (15) 

vl \MqJ VlOOkms^V 

2.3. Central wind 

We consider outflow rates < 10~^ Mq yr~'^ that can be driven by the stellar wind (Decampli 1981), as well as higher 
outflow rates < 10~^ yr"-'^ that can be driven by the interaction of the rotating magnetic fleld with the accreting 
disk (Shu et al. 2000; Konigl & Pudritz 2000). We do not consider the initial acceleration of the wind since it is 
expected to strike the disk surface after achieving terminal velocity. For simplicity, we assume a spherically symmetric 
wind with an isotropic distribution of density p^^ and velocity In this case, in a reference frame centered on the 
star, the wind density can be written as 

^-^W^^^^^ S^"^ 1^10-8 Mo yr-iji200 km s- J IrAuJ ' ^'^^ 

where = R'^ + and is the outflow mass loss rate per unit time. 

Our model could easily be extended to also include coUimated winds. We note that the assumption of spherical 
symmetry and isotropic distribution is required only in the small solid angle subtended by the disk. If the wind is 
substantially coUimated, then the true wind mass loss rate is substantially higher than the used in our "spherically 
symmetric" models. 
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Fig. 3. — Contours of the gas temperature (in units of 10"^ K) in the mixing layer surface density-disk radius parameter space. This is 
the surface temperature of the disk just below the mixing layer. 
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Fig. 4. — Mixing layer surface density as a function of disk radius for different boundary conditions and an entrainment efficiency e = 0.1. 
(a) Solutions for a wind incidence angle of 0.001° and different mixing layer surface densities at the inner boundary (lO"*^, 10"'', and 10"® 
g cm"'^, as labeled on each line), (b) Solutions for a mixing layer surface density of 10~* g cm~^ and different wind incidence angles at 
the inner boundary (0.001, 0.01, and 0.1°, as labeled on each line). 

3. RESULTS 

We find steady state solutions (i.e. 9tSm = dt^mVm = dt^m^m = 0) of equations (3)-(5) with the boundary 
conditions 5]m(-Rmm) — ^m{Rmin) — and j{Rmin, Zm{Rmin)) = 7oi where 7o is a small initial wind incidence angle 
at Rmin- We assume a solar mass star, and a ballistic wind with v^^ = 200 km s^^ and M^, = 10~^ M0 yr~^ in our 
fiducial model. 

We use the detailed thermo-chemical disk model of Gorti & Hollenbach (2008) to separately calculate the gas tem- 
perature, the sound speed, and the gas density at the boundary between ambient disk and mixing layer. This model 
includes FUV and X-ray radiation, dust collisions, photo-reactions, and chemistry. For simplicity, we extrapolate 
calculations at specific disk radii (1, 3, 10, 30, and 100 AU) to calculate the gas temperature at different radii. Figure 
3 shows contours of temperature in the disk radius-surface density parameter space. Here, surface density is measured 
from high z downward toward the disk midplane. The sharp transition at surface densities ~ 10~^ g cm~^ is due to 
dust extinction of stellar photons, and the transition at ~ 10~^ g cm~^ is caused by a complicated interplay of X-ray 
heating and ionization combined with [Ne II] and [Ar II] cooling. The transition at ~ 20 AU is due to a relatively 
sudden inability of FUV and X-ray heating to maintain gas temperatures ^ 300 K. Beyond this point [O I] 63 /im 
cooling can maintain gas temperatures < 200K. 

Figure 4 shows the mixing layer surface density for e = 0.1, a likely high efficiency, and different initial wind incidence 
angles and mixing layer surface densities at the inner boundary varying over several orders of magnitude. This figure 
illustrates that solutions at large radii (i? <^ 1 AU) are insensitive to the exact inner boundary conditions used for a 
given value of e. 

We summarize the main results for e = 0.1 in Figure 5. Although the wind-disk mixing layer accumulates material 
from the wind and the underlying disk on the way out, surface density variations in the mixing layer become small (at 
i? > 1 AU, see Figure 5a) due to the increasing surface area. 

The wind-disk mixing layer is above the disk scale height (Figure 5b) because the small wind incidence angles (Figure 
5c) result in correspondingly small wind ram pressures that can be balanced by the small disk pressure at large height. 
The disk pressure, Pd = prfC^(l -I- e^), is larger than the normal wind pressure, = PwV^ sin 7, (Figure 5d) and thus 
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Fig. 5. — Mixing layer and disk parameters as a function of disk radius for an entrainment efficiency e = 0.1. The inner boundary 
conditions are a wind incidence angle of 0.001° and a zero mixing layer surface density, (a) Disk surface density, E^, and mixing layer 
surface density, Em, as labeled, (b) Mixing layer height (solid line) and disk scale height (dotted line), (c) Wind incidence angle, 7, and 
mixing layer inclination, /3, as labeled, (d) Normal components of the disk pressure (solid line), = pdCs(l + e^), the wind pressure 
(dashed line), P„ = p„i;^sin'^7, and the curvature pressure term (dotted line), = Ti^v^ cos f^d/S/ dR. See equation (2) and Figure 2 
for a description of this pressure terms, (e) Mixing layer velocity (solid line), Vm, and escape velocity (dotted line), (f) Disk dispersal time 
(solid line), E^/S^a, and viscous evolution time scales (dotted lines) for different a viscosity values, as labeled. 

the mixing layer surface flares outward with increasing radius (Figure 5c). 

The mixing layer starts with roughly the same velocity as the wind at the inner boundary (Figure 5e) since the wind 
incidence angle is very small and very little disk mass has been entrained (Figure 5c); and slows down at larger radii 
as it accumulates mass from the disk (Figure 5e). Although the mixing layer velocity decreases with radius, it remains 
significantly higher than the escape velocity from the star. 

We compare the wind dispersal time (eq. [9]) with the viscous evolution time scale (eq. [13]) for a range of accretion 



rates in Figure 5f. For our fiducial model (M^, = 10 Mq yr 



200 km s ^, and e = 0.1) and an accretion rate 



'--^ 3 X 10~ M0 yr~ , viscous evolution is the dominant mechanism in the inner disk (< 50 AU) and wind stripping is 
the dominant mechanism in the outer disk in the absence of photoevaporation. We will explore the conditions under 
which wind stripping dominates the disk evolution in more detail below. 
We derive an analytic expression for the dispersal time by assuming that the curvature term contribution to the 
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Fig. 6. — Wind incidence angle (a), 7, and shock velocity (b), Vg = sin7, for an entrainment efficiency e = 0.1, stellar wind velocities 
in the range 100-300 km s~^, and wind outflow rates in the range lO^'^ — 10^** Mq yr~^, as labeled. Changing the entrainment efficiency 
results in negligible variations. 

normal pressure balance (eq. [2]) is similar to that of the stellar wind (Figure 5d). In this case, we can approximate 



the normal disk pressure as pdcj{l + e ) 



2pyjV^ sin^ 7, and the dispersal time (eq. 
TrEdCos /3csi?2(l + e2) 



as 



t 



where we use equation (IG) and r 
,(i?)~6Myr(^)"'(l + e2) 



MujW^esin 7 

R. We replace /? ^ 10° and 7^2° (Figure 5c) in equation (17) to obtain 
So \ f R\ f c, \ ( M, 
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g cm 
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AU 



4 km s 



10-8 Mgyr- 



200 km s' 



(17) 
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Note that the sound speed decreases with radius and mixing layer surface density (Figure 3). The decreasing sound 
speed means that Rcs does not increase much with R, which leads to a fairly constant dispersal time as a function of R in 
this case where surface density drops with R~^. Eq. [18] is in agreement with the estimate of HoUenbach et al. (2000). 
The highest dispersal rates (shortest dispersal times) correspond to the highest wind outflow rates and velocities, as 
expected. We can quantify the dispersing power of the wind with the parameter 



10-8 Mq yr- 



, 200 km s" 



(19) 



where we arbitrarily choose 77 = 1 for our fiducial parameters. 

Figures 6-8 compare results for models with different entrainment efficiencies, wind outflow rates, and wind velocities 
with our flducial model (e = 0.1, = 10-^ Mq yr-^, ti^, = 200 km s-^). These results agree well with the analytic 
approximation given in Eq. (18). 

We compare the wind incidence angles and the corresponding shock velocities for e = 0.1, wind outflow rates in the 
range 10~^ — 10"^ Mq yr~^, and wind velocities in the range 100 — 300 km s~^ in Figure 6. Variations due to changes 
in the entrainment efficiency are not significant. The wind incidence angle is not sensitive to e, M^, or Vw since the 
mixing layer is located several disk scale heights above the midplane, where large changes in pressure (and mixing 
layer inclination) can be accommodated by small changes in the mixing layer height. The jogs at ~ 20 AU are due to 
the rapid drop in gas temperature beyond this radius (Figure 3). We predict small wind incidence angles 7^4° and 
correspondingly small shock velocities (e.g. w^, sin7 — 17 km s-^ at 100 AU for AI^ = 10"* Mq yr-^ and — 300 
km s~^). 

We consider variations of the entrainment efficiency alone (e = 0.01 — 1) in Figure 7. Only the mixing layer velocity 
and the dispersal time are shown since variations of the other quantities shown in Figure 5 are not significant. Increasing 
the entrainment efficiency decreases the mixing layer velocity and reduces the disk dispersal time, as expected since the 
mixing layer receives more mass from the underlying disk in this case. Although the mixing layer velocity decreases 
with radius as it entrains more mass from the disk, it is possible to maintain the maximum possible entrainment 
efficiency, e = 1, since the mixing layer velocity remains larger than the escape velocity at all radii. 

Finally, we compare the dispersal mass flux rate and dispersal time due to wind stripping with those due to photo- 
evaporation and the mass flux and characteristic time due to viscous evolution in Figure 8. Since the dispersal time 
depends on the particular disk surface density distribution assumed, the mass flux outflow rate ("dispersal rate") is 
physically more meaningful. HoUenbach & Gorti (2009) show that EUV photons cannot penetrate the wind when the 
outflow rate is higher than ~ 10"^ Mq yr"^, while FUV and X-ray photons begin to penetrate once the outflow rate 
falls below ~ lO-'^ M© yr'^ Thus, we can ignore EUV photoevaporation for the high outflow rates (> 10 ^ Mq yr ^) 
considered here. We compare the dispersal rates by wind stripping with those due to photoevaporation driven by FUV 
and X-ray heating (Gorti & HoUenbach 2009, Figure 4). We assume outflow rate/accretion rate ratios x = M^/Macc 
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Fig. 7. — Mixing layer velocity (a) and disk dispersal time (b) as a function of disk radius for entrainment efficiencies in the range 
0.01 < e < 1, as labeled. The other wind parameters are standard (M^j = 10~* Mq yr~^ and = 200 km s~^). The dotted lines in 
panel (b) give viscous evolution time scales for different a viscosity parameters in the range 10--* <a < 10-2, as labeled. 

in the range ^ 0.01 — 1, as inferred from observations (White & HiUenbrand 2004, Figure 16, excluding edge-on disks). 
It is useful to consider different accretion rate regimes (Figure 8): 

• Low accretion rates (< 10~^ Mq yr~^. Figures 8a, b, c, and d). Viscous accretion is the dominant mechanism 
at most disk radii (< 150AU), while photoevaporation is the dominant mechanism in the outermost regions 
(> 150 AU). Wind stripping can be dominant at ^ 100 AU if the wind outflow rate approaches the accretion 

rate (x ~ 1). Once again, for Af^, < 10~^ M© yr~^, EUV photons can penetrate the wind (HoUenbach & Gorti 
2009) and EUV photoevaporation becomes dominant at ~ 1 AU (Liffman 2003; Alexander et al. 2006a). 

• Intermediate accretion rates (~ 10~^ Mq yr~^. Figures 8e and f). Viscous accretion is the dominant mechanism 
in the inner disk (i? < 60 AU) , while photoevaporation is the dominant mechanism in the outer disk. Since the 
accretion rate also depends on the surface density (eq. [14]), we assume a surface density higher than our fiducial 
value (Sd = 5 X 10^ g cm"^ at 1 AU) for intermediate and high accretion rates. 

• High accretion rates (> 10^® Mq yr^^. Figures 8g and h). Viscous accretion is the dominant mechanism at all 
disk radii. In this case, although optical photons can reach and heat the disk, FUV and X-ray photons cannot 
penetrate the correspondingly high outflow rates (HoUenbach & Gorti 2009) . 

4. SUMMARY AND CONCLUSIONS 

We present a model for the dispersal of protoplanetary disks by winds from either the central star or the inner 
disk. These winds obliquely strike the flaring disk surface and strip away disk material by entraining it in an outward 
radial- moving flow at the wind-disk interface, located at the surface of normal pressure balance several disk scale 
heights above the mid-plane. We derive conservation equations for the mass, momentum, and angular momentum in 
the wind-disk interface. 

The disk dispersal time scale depends on the velocity at which disk material is entrained into the mixing layer, which 
we quantify with an entrainment velocity ecs, where Cg is the sound speed and e is an entrainment efficiency. The 
dispersal time decreases as (ev^jM^)"^ (see eq. [18]). For disk surface densities dropping as R^^, the dispersal time 
is relatively independent of R (see Figs. 7 and 8 and Eq. [18]). We define a dimensionless parameter 77 cx ev^M^ 
that is unity for our standard parameters e = 0.1, Vyj = 200 km s^^, and M„ — 10~* Mq yr"^. The dispersal time 
for 77 1 (eq. [19]), tws{^ AU) '-^ 6(Eo/10'^ g cm"^) Myr, is somewhat larger than typical protoplanetary (dust) disk 
evolutionary time scales (~ 2 Myr) inferred from infrared observations (e.g. Haisch et al. 2001; Cieza et al. 2007). In 
fact, as discussed below, viscous evolution and photoevaporation hkely dominate disk dispersal. 

Figure 9 depicts the dominant physical mechanisms for different accretion rate regimes. For low accretion rates 
(< 10~* M0 yr~^), viscous accretion is the dominant mechanism at most disk radii (< 150AU), while FUV and X-ray 
photoevaporation is the dominant mechanism in the outermost regions (> 150 AU). In this case, wind stripping can 
be a dominant mechanism at ~ 100 AU only if the wind outflow rates approach the accretion rates, which is rare. 
EUV photoevaporation becomes a dominant mechanism at ~ 1 AU for low wind outflow rates (< 10~^ M0 yr"^). For 
intermediate accretion rates 10~^ Mq yr"^), viscous accretion is the dominant mechanism in the inner disk (< 60 
AU), while FUV and X-ray photoevaporation is the dominant mechanism in the outer disk. For high accretion rates 
(> 10~^ M0 yr"^), viscous accretion dominates in the entire disk. 

If wind stripping dominates viscous evolution, a gap may form if the inward accretion flow is less than the mass 
flux produced by wind entrainment at that radius, similar to the combined effects of viscous accretion and EUV 
photoevaporation (Clarke et al. 2001; Matsuyama et al. 2003; Alexander et al. 2006b). If the wind mass loss rate 
scales with the accretion rate, such as in X-wind models, then once a gap forms, accretion will rapidly decrease, the 
wind will turn off, and the gap will be refilled by viscous evolution of the outer disk. Time-dependent models which 
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Fig. 8. — Dispersal rate and the corresponding dispersal time (solid lines) for e = 0.1, = 200 km s ^, and different outflow mass loss 
rate/accretion rate ratios x = / Mace, as labeled. Dotted lines show the mass transport rate and the corresponding viscous evolution 
time scale due to accretion for accretion rates in the range 3 X 10~® — 10~® Mq yr~^ (and the corresponding a viscosity values using 
Eq. [14]), as labeled. Dashed lines show the photoevaporation dispersal rate and the corresponding dispersal time (from Figure 4 of 
Gorti & HoUenbach 2009) for the FUV luminosities, Lpyv, associated with the different accretion rates (Gorti &: HoUenbach 2008), as 
labeled. We assume a surface density higher than our fiducial value (S^ = 5 X 10^ g cm~^ at 1 AU) for correspondingly high accretion 
rates (e, f, g, and h). 
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(b) Low accretion rates (< 10 ^ Mg yi- ^) 



EUV Accretion Wind FUV 

— I 1 1 < 

1AU 10 AU 100 AU 



(c) Intermediate accretion rates (~ 10 ^ Mq yi- ^) 



Accretion FUV 
— I 1 1 

1AU 10AU 100AU 



(d) High accretion rates (> 10 Mq yi-"^) 



Accretion 

— I 1 1 

1AU 10AU 100AU 



Fig. 9. — Schematic illustrations of the physical processes that dominate disk dispersal (a); and the accretion rate and disk radii regimes in 
which viscous accretion, EUV photoevaporation (EUV), FUV and X-ray photoevaporation (FUV), and wind stripping (wind) are dominant 
mechanisms (b, c, d). The central wind interacts directly with the disk at all radii since the mixing layer flares outward (Fig. .5c). In the 
mixing layer at a given radius, disk material mixes not only with the wind but also with the wind material already mixed with the disk 
gas at smaller radii. Note that in (b), EUV dominates at ~ 1 AU only when the wind mass loss rate < 10~^ yr~^. 

include both the viscous evolution and mass loss by wind stripping are needed to explore this scenario. Once again, 
this requires wind outflow rates similar to the accretion rates (Figures 8a, c). 

Hartmann & Raymond (1989) found that shock heating at the wind-disk interface can explain the observed forbidden 
line emission from large radii (^ 50—100 AU). However, they predict shock velocities ~ 20 — 30 km s~^ (see their Figure 
6), while we predict smaller shock velocities (^17 km s~^. Figure 6b) and correspondingly small wind incidence angles 
(< 4°, Figure 6a), even for wind velocities of 300 km s^^, higher than the maximum value (200 km s~^) assumed by 
Hartmann & Raymond (1989). The difference arises because Hartmann & Raymond (1989) assume that the tangential 
momentum of the mixing layer is given by the total accumulated tangential momentum of the wind while we solve for 
tangential momentum conservation [eq. (4)]. Note that mass input from the disk affects the tangential momentum of 
the mixing layer since it increases the mixing layer surface density. The observed [Ne II] 12.8 /im forbidden line emission 
(e.g. Lahuis et al. 2007) may be explained by EUV or X-ray heating (Glassgold et al. 2007; HoUenbach & Gorti 2009), 
although the high [O I] 6300 A luminosities observed in some cases (Hartigan et al. 1995) defy explanation. The 
combination of EUV heating and ionization combined with turbulent mixing and heating in the wind-induced shear 
layer may provide the high gas temperature (> 3000 K), high electron abundances, and high atomic O abundances 
needed to explain the [O I] luminosities. Further work is needed in this area. 

We thank Uma Gorti for supplying us with prepublication results of her disk models that allowed us to estimate 
the temperature structure of the disk just below the shear layer created by the wind. This research was supported 
by the Miller Institute for Basic Research in Science, University of California, Berkeley; the Carnegie Institution of 
Washington; a Natural Sciences and Engineering Research Council of Canada grant; the NASA Astrophysical Theory 
program; and the NASA Astrobiology Institute. 
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APPENDIX 
CONSERVATION EQUATIONS 

We consider normal pressure balance; and mass, momentum, and angular momentum conservation for a mixing layer 
annulus between R and R+ SR, as shown in Figure 1. We denote the wind density, velocity, and incidence angle as 
Pw, Vw, and 7 respectively; the disk density and sound speed as pd and Cs respectively, and the mixing layer surface 
density, velocity, and inclination as S^, Vm, and (3 respectively. 

Force balance in the direction perpendicular to the surface of the annulus can be written as 

2tt R5spy,vl^ siv? 7 + 2nRj:„,vf^ sin 6(3 = 2TrR6spdcl + 2nRSspde'^cl. (Al) 

On the left-hand-side (LHS) of equation (Al), the first and second terms correspond to the force associated with the 
normal momentum flux of the wind and the normal momentum flux through the inner boundary respectively (Figure 
1). On the RHS of equation (Al), the first and second terms correspond to the thermal and ram pressure forces of 
the disk respectively. Replacing Ss — SR/ cos/? and taking the hmits SR^ and 5/3 ^ in equation (Al) yields 

dd 

p^vism^j + i:^vf^cos(3—=pdcl{l + e^). (A2) 

The second term on the LHS of equation (A2) is given by the product of the tangential momentum flux per unit length 
and time, Em^^mi ^^'^ the curvature of the mixing layer surface, cos PdP/dR = dp/ds. 

The mass and momentum in the direction tangent to the surface are 2TTRSsT,m and 2TrRSsT,mVm respectively. The 
mass change per unit time is given by mass input from the wind and the disk, and mass transport through the annulus: 

d 

— (27ri?S„i(5s) = 2ttR5sPwVw sin 7 + 2TrRSspdeCs + 27ri?S„iW,„ — 27r_R^S^jW^j, (A-3) 

where we define i?+ = R + 6R, S+ = E„j(i?+, z(i?+), t), and w+ = Um(i?+, z{R'^), t). The last two terms on the 
RHS of equation (A3) correspond to the mass flux per unit time across the inner and outer boundaries of the annulus 
respectively. Replacing 5s — 6R/ cos/3 (Figure 1) and taking the limit SR — > in equation (A3) yields 

. cos /3 9 , , , , ^ , 

—g^ = PwVyj Sm 7 -I- PdeCs — ^ {RY.mVm ) ■ (A4) 

Momentum conservation can be written as 

— (2TTRT,,nSsv,n) = 2-KR5spwv'^ sin 7 COS 7 -I- 2'KRY,mv'L cos 6(3 — 2'KR^T,t,vt^ , (A5) 
ot 

where the last two terms on the RHS are the momentum flux per unit time across the inner and outer boundary of 
the annulus respectively. Replacing Ss = SR/ cos/3 and taking the limits S/3 and SR ^ in equation (A5) yields 

d cos /3 d 

— {T,„iVm) = Pwvl sin 7 cos 7 — — (i?S„t)^,) . (A6) 

Assuming that the wind does not carry significant angular momentum, angular momentum conservation can be 
written as 



— {2TTR^,nSsn„,R^) ^2nRSspdeCsnkR^ + 2^RY.mVrr.^mR' - 2^i?+S+w+r!+i?+2, (A7) 

where fi^ = {GM / R^Y^"^ is the Keplerian angular velocity of the disk and ri„i is the angular velocity of the mixing layer. 
The angular velocity of the disk must be sub-Keplerian due to radial pressure support, by a factor ~ 1 — 0{H/R)^, 
where H is the disk scale height. We ignore this effect since it is negligible for protoplanetary disks with (H/R)'^ <C 1. 
The angular momentum loss due to the outward motion of the mixing layer [third and fourth terms on the RHS of 
equation (A7)] is partially compensated by the angular momentum input from the underlying disk [first term on the 
RHS of equation A7]. Once again, replacing Ss — SR/ cos (3 and taking the limits Sf3 ^ and SR ^ in equation 
(AT) yields 

d cos /3 d 

— (S^rim) = PdeCsfl H3~^ {R^T,jnVm^jn) ■ (A8) 
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